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THE AXIALLY SYMMETRIC RESPONSE OF AN ELASTIC CYLINDRICAL 


SHELL PARTIALLY FILLED WITH LIQUID 

By Richard M. Beam and LeRoy R. Guist 
Ames Research Center 


SUMMARY 


The axially symmetric response of a circular cylindrical shell partially 
filled with a liquid and subjected to axial displacements and pressurization 
fluctuations is investigated analytically and experimentally. Two solutions 
are presented: one, an "exact 11 series solution and the other, an approximate 

solution. The correlation of theoretical and experimental results indicates 
that the theory presents a good mathematical model of the physical system and 
that the approximate solution is adequate for many practical problems. 


INTRODUCTION 


The use of liquid propellant for missiles and spacecraft launch vehicles 
introduces problems associated with the dynamics of liquid-filled shells. 
Initially, there were primarily control problems resulting from the "propel- 
lant sloshing" or motion of the liquid free surface. Later, however, as the 
size of propellant tanks was increased and the ratio of wall thickness to 
radius decreased, the coupling of the shell and liquid dynamics produced 
unexpected and undesirable oscillations of flight vehicles. In particular, 
the longitudinal oscillation resulting from the coupling of the liquid-shell 
system with the shell support structure and the tank pressurization system 
can influence the success of certain missions and can Interfere with the 
comfort of astronauts in manned vehicles . 

Several papers have been published on the dynamics of the coupled liquid- 
shell system (refs. 1-8). Most of these papers are limited to the anaxisym- 
metric oscillations generally predominant in flutter calculations. Those 
papers that consider the axisymmetric oscillations lack information required 
to analyze the complete vehicle, namely, (a) the equations and coefficients 
that couple the liquid-shell system and the shell support and pressurization 
systems, (b) approximate equations suitable for preliminary design calcula- 
tions, and (c) experimental verification of the applicability of the theory. 

The purpose of this report is to present a method for computing the 
response of a preloaded, thin circular cylindrical shell partially filled with 
a liquid and subjected to axial base displacements and applied pressure fluc- 
tuations in the liquid-free shell space. A general series solution is 
obtained for the natural frequencies of oscillation and normal modes . In 
addition, an approximate one -term solution applicable to very thin shells is 
presented. The equations coupling the motion of the liquid-shell system with 



the base and pressure fluctuations are derived. Finally, a comparison is 
made between the two analytical solutions and experimental results. 
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SYMBOLS 

constant defined by equation (19) 
coefficient of series of modal functions 
constants defined by equation (19) 
constant defined by equation (29c) 
constant defined by equation (19) 
modulus of elasticity of shell material 
dimensionless length of shell, 

R o 

depth of liquid in shell 

modified Bessel function of the first kind 

Bessel function of the first kind 

integers 

length of shell 

integers 

shell axial preload (i.e., axial shell loading for static 
condition) 

N x 

dimensionless shell axial preload, 

force on bottom of shell 
constant defined by equation (19) 

pressure in region occupied by liquid (also used as integer 
subscript) 

pressure in liquid-free volume of shell 

Po 

dimensionless pressure in liquid-free volume of shell, g— 
constant defined by equation (l 9 ) 
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z in^2n; z 3n 

a 
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Pk 
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5 n 
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T 1 
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surface pressure on shell 
cylindrical coordinate 

R 

dimensionless cylindrical coordinate, — 

K o 

radius of cylindrical shell 
constant defined by equation (19) 
coefficient of Fourier -Bessel series 
integer 

coefficients of Fourier sine series 
time variable 

axial displacement of shell bottom 

dimensionless axial displacement of shell bottom 

radial displacement component of shell wall 

dimensionless radial displacement component of shell wall, 
w 

Ho 

cylindrical coordinate 

coupling coefficients defined by equations (29a), (29b), and 

(33) 

x 

dimensionless cylindrical coordinate, — 

R o 

dimensionless shell thickness, 

K o 

constant defined by equation ( 6 b) 
dimensionless form of 5 n , s n Ro 
shell thickness 

constant defined by equation ( 6 a) 
dimensionless liquid depth, 

R o 

integration variable 

step function defined by equation (l5a) 
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P* 

v 

I 

P 

Ps 

or 

cT 

0 

9s 

? 


Ps 

ratio of shell -material and liquid densities, — 

Poisson 1 s ratio for shell material 
integration variable 
liquid mass density- 
shell- material mass density 
eigenvalue (frequency squared) 
dimensionless form of a, a ■ ° P 

velocity potential function 
modal function 

displacement potential function 
frequency of oscillation 

AmLYSIS 


The liquid-shell geometry and coordinate system 'are shown in figure 1. 

The analysis is restricted to small (linear) displacements that are symmetric 
about the axis of the cylinder. The liquid is incompressible and inviscid 
and the flow irrotational. The liquid is bounded at the bottom of the cylin- 
der by a rigid flat surface, at the free surface by a light (negligible den- 
sity) gas at pressure p G (t), and at the cylinder surface by the flexible 
shell wall. The shell has simply supported boundaries and is preloaded with 
a constant axial stress. Only the radial inertia of the shell is included in 
the analysis. All body forces acting on the liquid are assumed to be zero; 
therefore, the "sloshing” modes of the free surface do not appear in the math- 
ematical model. In general, the body forces may be excluded from this type 
of analysis because of the low frequency spectrum of the sloshing modes 
relative to the liquid shell modes. 

The analysis is divided into four parts: (a) equations of motion of the 

liquid, (b) coupling of the shell -liquid equations, (c) determination of natu- 
ral frequencies and modes of oscillation, and (d) response of the liquid-shell 
system to displacement and pressure disturbances. 


Equations of Motion of the Liquid 

For irrotational flow of an incompressible liquid, the velocity potential 
0 satisfies the continuity equation (ref. 9) 
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0 = 0 


( 1 ) 


v 2 <d 



1 , 
R 3R + 



in the region occupied by the liquid (0<x<h, 0 < R < Ro) • Since the 
liquid displacements at the shell are directly related to the shell wall dis- 
placements , it is convenient to use a displacement potential, Y (where 
$ = S^/St), rather than the velocity potential. In addition to being a har- 
monic function (v 2 ^ =■ 0) , Y must satisfy the following boundary conditions. 

On the free surface (x = h) , the pressure is a prescribed function 
[p Q (t) ] time 


P 


x=h 



x=h 


= p 0 ( t ) 


( 2 ) 


where p is the pressure in the liquid and p is the liquid density. 

At the shell bottom (x = 0) , the liquid must have the same axial motion 
[u(t) ] as the rigid tank bottom 


By 

Sx 


= u(t) 

x=o 


( 3 ) 


The liquid adjacent to the shell wall (R = R 0 ) must move radially with 
the tank wall 


SY 

Sr 


R=Ro 


w(x,t) 


( 4 ) 


where w(x,t) is the shell -wall radial displacement component. 

The displacement potential Y may be obtained by the method of separa- 
tion of variables applied to the continuity equation (v 2 ^ = 0) . If the arbi- 
trary constants resulting from the separation of variables are chosen so that 
boundary condition equations (3) and ( k ) are satisfied, the following function 
is found 
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¥(x,R,t) = r Q (t) + V r n (t) ~y \ + -“r- ( r2 - 2X 2 w(^t)d| 

° ^ n 5 n smh B n h J 0 (6 n R 0 ) 2hR 0 J o 

n=i 

0° 

- 1 (f w<6 ’ t)cos ^ 4 0 k- »« ^ + “(*) <5) 

k=i 


where 


Jo( 8 n R o) =0 n = 0, 1, 2 


Pk = I 1 k = 1, 2 


(6a) 

(6b) 


and J Q and I Q are Bessel functions and modified Bessel functions, respec- 
tively, of the first kind and zero order. The arbitrary functions of time, 
r n , may be used to satisfy the remaining boundary condition (eq. (2)). 

The value of r n may be determined in the following manner. The func- 
tion ¥ is substituted from equation (5) into the free surface boundary con- 
dition equation (2). The resulting equation is then multiplied by RJ 0 (s s R)dR 
and integrated from zero to R Q . After the appropriate algebraic manipula- 
tions, the following expressions for r n result 1 


1 The derivation of equations (7) and (8) is expedited by the following 
identities 

(a) cos P^h = | x 2 cos p k x dx 
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J w(£,t)cosh S s £ d£ 


( 7 ) 

(8) 


The pressure in the fluid may now be expressed as a function of the 
shell-wall displacement , the free-surface pressure, and the axial motion of 
the shell bottom: 
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( 9 ) 


It is worthwhile to note that the expression for the pressure does not 
depend on the shell boundary conditions and, accordingly, is applicable for 
the analysis that includes any shell boundary conditions. 

The liquid pressure, q, acting on the shell is obtained from 
q = p(x,R 0j t) 
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With the dimensionless parameters. 
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equation (l0) may be rewritten 
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Equations of Motion of the Shell 


If a circular cylindrical shell of uniform thickness, 5, is subjected to 
a constant axial load , 2 * N x , uniformly distributed around the circumference, 
and an axially symmetric surface load, q(x,t) , the differential equation for 
the radial deflection component w (ref. 10 ) is 
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12(1 - v 2 ) 3x 4 


d 2 w , E5 aw e ^ 

- N x ~s. + r-2 w = -P s r-2 5 + 

OX h 0 ot 


(13) 


where E, v, 5, and p s are the modulus of elasticity, Poisson 1 s ratio, 
thickness, and density of the shell material. 


With the dimensionless notation 



_ N. 
N x = 
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and q(x,t) from equation ( 12 ) , the shell equation may be written 
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(15) 


2 The value N x is determined from the static loading condition, that 

is, from any constant pressure in the shell plus axial loads applied to the 
ends of the shell. It is assumed that the variation of axial loading result- 
ing from p Q (t) and u(t) produces negligible variations of N x (i.e., 

N x - constant) . 
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where 


A (a) 


1 0 < a < e 

0 a > G 


(15a) 


For simply supported boundary conditions, the shell displacement and its 
derivatives must satisfy the following equations: 


w(0) = w"(0) = w(e) - w TT (e) - 0 


(16) 


where 




w 


w 1 


Shell -Liquid Natural Frequencies and Modes of Oscillation 

The shell -wall response to base displacements and pressure fluctuations 
is the function w which satisfies equations (l5) and (l6) plus the initial 
conditions. This response, which can be used to evaluate the desired cou- 
pling characteristics is obtained in the following manner. First, a series 
solution for the homogeneous part of equation (15) (p Q = u = 0) will be 
obtained by assuming that the time variation of w is proportional to 
exp( t) . This will result in a set of eigenfunctions cp s (<x) that can be 

used in computing the response w(a,t) to the disturbances p Q and u (i.e., 
the nonhomogeneous eq. (15))- In this manner, let 
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00 

^ sj sin a exp(i \fo t) 


(17) 


and substitute into equation (15) • The equations for the coefficients sj 
are obtained by multiplying the resulting equation by sin(mjt/e)a da and 
integrating from zero to e. The following set of simultaneous equations 
results : 
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With the definitions 
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equation (l8) can be simplified to 
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The simultaneous equations ( 20 ) are homogeneous and therefore have a non- 
trivial solution only if the determinant of the coefficients of sj vanish. 
This provides the eigenvalue problem for the eigenvalues a s and the cor- 
responding eigenfunctions sj s which, in turn, lead to the modal functions 


^s 


00 

Sj s sin a 


( 21 ) 


by 


The natural frequency of oscillation of the mode cp s is related to a s 


to 


s 



( 22 ) 


where to g is in radians per unit time. 


Response to Pressure and Displacement Fluctuations 


With the aid of the 
the solution to equation 


modal functions obtained in the previous section, 
(15) may be expressed as 


w(x,t) 
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a n (t)cp n (a) 

n=i 


(23) 


where a n is as yet an undetermined function of time. The differential 
equations for evaluating a~ are obtained by substituting w from equa- 
tion (23) into equation ( 15 ), multiplying the resulting equation by cp m (a)da, 
and integrating from zero to e, which produces 
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With the aid of equations (15) and (l6) and the properties of the eigen- 
functions of the previous section, the following orthogonality relationships 
can be derived: 
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The differential equation for a n (t) (eq. (24)) may be simplified to 


where 


a n (t) + w n e a n (t) = Z in u(t) + Z 2nPo (t) 
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If the disturbances u and p Q are prescribed, the shell response is 
readily obtained from equations (23) and (28) . The coupling of the liquid- 
shell system and the shell- support system requires, in addition, the axial 
force acting on the bottom of the shell. The pressure in the fluid as a func- 
tion of a n is obtained by the substitution of w from equation (23) into 
the pressure equation (9) • 
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The total force, P, acting on the bottom of the shell is J* p(o,R,t)dA, where 
the integral is over the area of the shell bottom; therefore, 
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The coupling of the liquid-shell system and the shell-support and pressuriza- 
tion system is thus provided by equations (28) and (32). The coefficients 
Z x , Z 2 , and Z3 will be called "coupling coefficients" in the discussion to 
follow. 
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Approxima te solution (thin shells) . - The complexity of the frequency 
equation and coupling coefficients in the previous sections makes that solu- 
tion unattractive for preliminary design or analysis of a system. Since most 
liquid-shell systems used in missile and spacecraft construction have very 
thin wall shells, an approximate membrane solution of respectable accuracy can 
be expected. For thin shells subjected to loading that changes only gradually 
over the shell surface, the bending rigidity may be neglected in the interior 
of the shell. In addition, for the liquid-shell system under consideration, 
the shell boundary conditions (eq. (l6)) may also be neglected since their 
influence is localized near the boundary. 


A displacement potential satisfying the fluid boundary conditions 
equations (2), (3), and (k) is 


-u ~ 2 w(£,t)cos d£ 
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where 
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^k ~ 27 h 


(6b* ) 


Note that as defined by equation (6b * ), differs from the 3^ of equa- 

tion (6b) . The simplification of the displacement potential and thus the new 
definition of (3 results from the relaxation of the boundary conditions on 
the shell and the lack of shell bending rigidity, that is, the displacement w 
at x = 0, is not required to be zero and, in addition, it is known a priori 
that the shell displacement, w, will be zero at the free surface (x = h) . 

The pressure in the liquid is 
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and the surface loading q(x,t) on the shell is 
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The equation for the shell- wall displacement w (neglecting shell bending 
rigidity) is 
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The frequency equation is obtained by substituting 
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into the homogeneous part of equation (15’) (i-e., p o =u=0) 
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The frequency equation is thus reduced to a single equation as opposed to the 
infinite determinant of the previous solution where the shell boundary 
conditions were included. 


The response of the shell wall can be expressed in terms of the modal 


functions cp . = cos 
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The equations for a n are obtained by substitution of w from equation (23 l ) 
into equation (l5 T )j multiplication of the resulting equation by 


cos 
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da and integration over the interval zero to e 
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where 
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The total force on the shell bottom P = j* 1 p( 0 ,R,t)R 6B dR is 
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The frequency eqioation and coupling coefficients of this approximate 
solution are obviously more tractable than those of the series solution. The 
accuracy and applicability of the two solutions will be discussed further in 
the succeeding sections. 


COMPARISON OF SERIES AND APPROXIMATE SOLUTIONS 


This section presents the results of a brief numerical investigation 
designed to illustrate the convergence of the series solution and to compare 
the series and approximate solutions. First, the convergence of the series 
solution with increasing size of the frequency determinant is illustrated and 
compared with the approximate solution. Next, the convergence of the series 
solution with increasing shell thickness -to -radius ratio is considered, and, 
finally, the mode shapes and coupling coefficients from the two solutions are 
compared. Shells with large radius -to -thickness ratios (i.e., thin shells) 
were selected for the examples since the series solution rate of convergence 
(with frequency determinant size) should be slowest in these instances. The 
following discussion is limited to the first (lowest) two frequencies since 
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these are generally of primary concern to the analyst. The accuracy of the 
two solutions (as compared with experiment) will be discussed in the next 
section. 


Frequency Parameters 

The series solution frequency equation (determinant of coefficients 
of eq. (20) equated to zero) was approximated with various numbers of rows 
(m) and columns (j) (m = j) for a liquid-shell system with the following 
properties: 


e = 2. 4 

e = 2.0 

N x = 0. 286x10 “ 6 
v = 0.3 

p = 8.6 

p = 0.00044, 0.044 

The eigenvalues and eigenvectors of the resulting frequency equation 
were obtained with the aid of a digital c omputer using the threshold Jacobi 
eigenvalue routine. The variation of frequency parameter with determinant 
size is shown graphically in figure 2. The first and second mode frequency 
parameters (cf/p) are shown for two values of shell thickness -to -radius ratio 
(P) . Note that the second mode frequency "jumps" as the size of the determi- 
nant is increased from M = 11 to m = 13 and, correspondingly, the perturba- 
tion of the first mode frequency parameter plot is slight (almost undetectable 
in the figure) in the same region. Examination of the mode shapes reveals 
that this discontinuity is possibly a result of convergence to an "incorrect" 
second mode (i.e., lack of sufficient terms in the series to adequately 
describe the second mode) if m < 13 . 

For the thinner shells (P = 0.00044) the gradient of the shell- wall dis- 
placement should be very large in the vicinity of the shell bottom and, there- 
fore, a slow convergence of the mode shape and frequency would be expected 
(due to the behavior of the sine series solution for a = 0) . As the shell 
thickness is increased, the convergence should improve. This is borne out 
by the frequency parameter plot (fig. 2) for the shell 100 times the thickness 
(p = 0.04) of the thinner shell. The approximate solution for the two shells 
is shown in figure 3* Note that for the thinner shell the difference between 
the two solutions is less than 2 percent for the first mode and 6 percent for 
the second mode. Even for the thicker shell (p = 0.04) the discrepancy 
between the series and approximate solution frequency parameter is less than 
5 percent for the first and second modes. 
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MODE SHAPES 


The mode shapes corresponding to the previous shell parameters are shown 
in figure 3- (The results for the approximate solution are indistinguishable 
for p = 0.044, p = 0.00044.) Here again the improvement of convergence of 
the series solution with increasing shell thickness is demonstrated. The 
first mode shape is, in fact, probably better represented by the approximate 
solution than by the series solution for a > 0. 


Coupling Coefficients 

The natural frequencies and coupling coefficients are or primary impor- 
tance in the analysis of the liquid-shell system. The coupling coefficients 
for the two shells under discussion are plotted in figure 4 as functions of 
frequency determinant size. (The results for the approximate solution are 
indistinguishable for 3 = 0.044, P = 0.00044.) Note the relatively small 
percent discrepancy (less than 6 percent except for Z 22 ) between the series 
and approximate solutions for the thinner shell and the only slightly 
increased discrepancy for the thicker shell. 

This rather limited numerical investigation indicates that the conver- 
gence of the series solution can be expected to be good even for thin shells 
and that the approximate and series solutions differ by small percentages for 
thin shells and are, therefore, almost equally applicable. As the shell 
thickness increases (3 increases), the discrepancy between the two solutions 
increases and the convergence and applicability of the series solution improve 
correspondingly . 


EXPERIMENTAL OBJECTIVES 


The experimental work was conducted to determine the dynamic properties 
of very thin liquid-filled cylindrical tanks subjected to axial vibration. 
The specific objectives were: 

1. To determine the fundamental natural frequency of axisymmetric 
motion and the natural frequency of as many higher modes as feasible for 
various liquid levels and internal pressures, 

2. To determine the frequency response curve for the fundamental axi- 
symmetric mode indicating the ratio of internal pressure at tank bottom to 
tank base acceleration versus driving frequency, 

3. To identify the asymmetric modes of vibration at several values of 
internal pressure and water level. This was done primarily to aid in dis- 
tinguishing the symmetric modes from the asymmetric modes. 
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Test Specimens 


The test specimens used in this study are shown in figure 5* The first 
specimen was a 5 -inch- diameter, 6-inch-long, 0 .001 -inch -thick brass shim stock 
tank which was later replaced by a 0 . 0007 -inch -thick electrodeposited nickel 
tank of the same dimensions. An 11 -l/2 -inch -diameter, 23-1/2 -inch -long, 

0 .001 -inch -thick brass shim stock tank was us.ed in the final phase of testing. 
The specimen size and thickness affected the choice of test apparatus and 
equipment. The following factors were given primary consideration in 
selecting specimen and apparatus properties. 

1. Very thin tanks are difficult to manufacture and handle, but have 
lower natural frequencies for a given tank size, or for a given natural 
frequency, the tank size is smaller. 

2. A large tank size requires a larger shaker and support fixture, 
resulting in lower frequency elastic modes within the shaker -fixture system. 

3. A large thin tank with the resulting lower natural frequency and 
higher mass permits better control of acceleration level and waveform at low 
acceleration levels. 

k. A large thin tank has a low burst pressure; hence, permissible 
acceleration levels are low. 

5. A low-frequency test range provides a better acceleration waveform 
from the shaker and less distortion of the waveform due to shaker -fixture 
dynamics . 

The 11 -l/2 -inch -diameter specimen was selected as the best trade off 
among the above factors. The smaller tanks used in earlier tests were easier 
to fabricate and were intended as preliminary models to indicate problem 
areas in testing technique and data quality. 

As shown in figure 5, the specimens were constructed with thick end 
flanges so that they could be attached to a base plate and fitted with a lid. 
The cylindrical shell was soldered directly to the end flanges, providing 
essentially a "built in" end condition, since it was not feasible or necessary 
to simulate the "pinned" end condition assumed in the analysis for such a 
thin shell. The base plates were designed to support the tank on the vibra- 
tion fixture, to provide a pressure seal, and to carry pressure cells. The 
lids also served as pressure seals, and means for filling and pressurizing 
the tank. 


Test Apparatus 

The general arrangement of the test apparatus is shown in figure 6. An 
electrodynamic shaker driving a hydraulic "slip table" produced the unidirec- 
tion sinusoidal motion of the specimen. The tank was pressurized with a shop 
air supply through a mercury manometer. Before the tank was pressurized, the 
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water level was measured with a water manometer connected to the tank base. 
The weight of the support fixture, the moving element of the slip table, and 
the specimen were counterbalanced by means of a hanging weight and pulley 
arrangement using elastic "bungee" cord to isolate fixture motion from the 
weight . 


Instrumentation 

The instrumentation arrangement is shown in figure 7* A differential 
pressure cell in the center of the tank base plate sensed the difference 
between pressures at the bottom of the liquid and the airspace above the 
water. Two accelerometers at the tank base provided signals for controlling 
shaker amplitude and for measuring base acceleration. The use of dual accel- 
erometers eliminated any interaction problem between the recording instruments 
and the shaker servo-control circuit. Lateral motion of the tank wall was 
detected by capacitive type noncontacting displacement probes. These probes 
were attached to an overhead mounting fixture that provided axial and rotary 
movement of the probes to any point on the tank surface. 


Procedure 

The test procedure was essentially the same for all specimens. First, 
the frequency of longitudinal motion was swept from 50 to 1200 cps for the 
5 -inch -diameter tanks and from 50 to 500 cps for the 11 -l/2- inch -diameter tank 
for a number of water levels. These sweeps were intended to locate and iden- 
tify as many system modes as possible. The acceleration level was varied ini- 
tially to determine, by trial and error, a level for which subharmonic surface 
resononce and liquid surface instabilities (ref. 11 ) produced minimum distor- 
tion of acceleration and pressure signals. The upper limit of the accelera- 
tion level was of course determined by the requirement that peak pressures 
encountered at symmetric mode resonance be below the tank yield pressure. 

The amplitudes of base acceleration, center pressure, and top accelera- 
tion, as well as phase angle relative to base acceleration, were recorded 
during each frequency sweep. Frequency sweeps were made for various values 
of internal pressure up to 1.5 psi for the large tank and k psi for the small 
tanks, and for various water levels from empty to full. Asymmetric modes 
were identified by rotating the displacement probes around the tank and noting 
the variation of wall motion amplitude as a function of angular position of 
the probe while frequency and base acceleration amplitude remained constant. 
These modes are easily identified since they characteristically exhibit an 
even number of peaks of radial displacement amplitude in 3^0° of rotation. 

First symmetric mode frequency response data were obtained by sweeping 
the excitation frequency at a constant acceleration and observing base accel- 
eration and pressure amplitudes. The response curve for the 5 -inch -diameter 
tank was taken for the nearly full condition (approximately l/2 inch below 
the tank lid), while the half full condition was used for the 11-1/2-inch- 
diameter tank. 
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RESUHJ?S 


The first symmetric mode frequency is plotted in figure 8 as a function 
of water level for 3*0 psig internal pressure for the 5 -inch -diameter nickel 
tank. Figure 9 is a similar plot for the 11 -l/2 -inch -diameter tank with 1.0 
psig internal pressure. In both figures, the theoretically predicted curves 
are shown for comparison . Figure 10 is a plot of the first symmetric mode 
frequency as a function of internal pressure for the 5 -inch -diameter tank 
nearly full (approximately l/2 inch of airspace above the water) . 

The experimental frequency response curves for the 5 -inch-diameter and 
11 -l/2 -inch -diameter tanks, respectively, are shown in figures 11 and 12. 

These curves show the ratio of pressure amplitudes at the center of the tank 
base to base acceleration as a function of driving frequency. For the 5 -inch 
tank, the nearly full condition with 3 *0 psig internal pressure was chosen, 
while a 12-inch water level and 1.0 psig internal pressure was used for the 
11 -l/2 -inch- diameter tank response curve. In either case, water level, inter- 
nal pressure, and acceleration level were chosen to yield the best data. 

Test conditions were selected by trial and error to obtain the set of condi- 
tions where acceleration and pressure waveforms showed the least distortion 
from subharmonic surface resonance, liquid surface instabilities, asymmetric 
breathing modes, and shaker -fixture dynamics. The extent of distortion 
appeared unpredictable, but certain combinations of test conditions yielded 
better waveforms than others. 

These distortions were clearly caused, in part, by the subharmonic 
liquid surface resonance and capillary wave formation associated with breath- 
ing modes as reported by Kana, Lindholm, and Abramson (ref. 11) . They also 
reported that, at certain acceleration levels, a jump phenomenon caused a 
sudden change in breathing mode response amplitude when the natural frequency 
was approached from below. The beat phenomenon (ref. ll) that involved a low 
frequency amplitude modulation of the pressure and radial wall displacement 
signals was also observed at certain frequencies. This amplitude modulation 
was always accompanied by a large, visible liquid surface motion and a modu- 
lation of the audible energy from the tank wall, both at the same frequency. 

No attempt was made to study these phenomena quantitatively because the test 
apparatus was inadequate for such an analysis and they were not important to 
the purpose of the present study. 

The frequency spectrum of the most predominant asymmetric breathing 
modes for the 5 -inch -diameter brass tank is plotted in figure 13- These 
modes were identified to determine when they were likely to distort the 
symmetric mode response. 


Discussion of Experimental Results 

The comparison of experimental and analytical frequency versus water 
level curves (fig. 8) indicates agreement within 10 percent for water levels 
above 2 inches. Below the 2-inch water level, however, the experimental 



values of fundamental frequency are as much as 1 6 percent higher than pre- 
dicted. For water levels above 3“l/2 inches , both fundamental and second 
mode data agree with predicted values within 2 percent. 

The comparison of predicted and experimental frequency versus water 
level curves for the 11 -l/2 -inch -diameter tank (fig. 9) indicates agreement 
within 5 percent for all data obtained. Data for water levels above 14 inches 
were not obtained because the test specimen ruptured accidentally at that 
point . Attempts to locate the fundamental frequency at water levels below 
k inches were unsuccessful because the distortion of both acceleration and 
pressure signals was high throughout the frequency spectrum. Testing with a 
new specimen at higher water levels was discontinued since the data of both 
figures 8 and 9 indicated very good agreement between prediction and experiment 
at higher wat er levels . 

The measured variation of fundamental natural frequency with internal 
pressure (fig. 10) indicates a frequency change of approximately 2 percent as 
the pressure varies from 1 to 4 psi, which agrees with the theoretical result 
that frequency is essentially independent of internal pressure (see eq.(20 ! )). 
Data from the other two specimens also indicate that frequency is essentially 
independent of internal pressure. 

The frequency response curves of figures 11 and 12 indicate that the test 
system behaves essentially as a linear single-degree-of -freedom system near 
the fundamental frequency. The amplification factor was not calculated since 
no attempt was made to assess system damping. 

The natural frequencies of symmetric modes of the liquid-tank system 
were difficult to determine primarily because of the presence of other phe- 
nomena (ref. 11) that distorted the symmetric mode motion. The excitation of 
asymmetric tank breathing modes at frequencies very near symmetric mode fre- 
quencies was almost unavoidable partly because these breathing modes are 
capable of subharmonic excitation. Large amplitude breathing mode motion at 
l/2 and l/3 harmonics was not uncommon and it is suspected that even lower 
harmonics were present although they could not be clearly identified. The 
fact that breathing mode frequencies vary with pressure, however, while the 
symmetric mode frequencies are essentially independent of pressure provided a 
means of isolating the symmetric modes. Another phenomenon that caused diffi- 
culty, primarily in attempting to establish precise pressure-to-acceleration- 
amplitude ratios, was the frequency modulation of symmetric mode motion at 
the same frequency as some low frequency liquid surface mode. This phenom- 
enon was more pronounced at higher acceleration levels and resulted in large 
surface motions. A difficulty which arose from the test apparatus configura- 
tion was that of producing a pure sinusoidal axial motion of the tank base as 
a rigid plate. This became impossible at higher test frequencies where the 
inevitable deformations of the base undoubtedly excited unwanted system modes. 
As mentioned earlier, this problem would be difficult to remedy because, even 
though larger tanks have lower frequencies, base flexibility also increases 
with tank size. 
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CONCLUDING REMARKS 


A comparison of the two theoretical solutions demonstrated that the 
series and one-term approximate solutions are equally applicable for thin 
shells; however, the ease of application of the approximate solution warrants 
its use. As the shell thickness -to-radius ratio becomes larger, the series 
solution converges more rapidly while the accuracy of the approximate 
solution deteriorates. 

The experimental results of this study indicate that the theory does, in 
fact, provide a good mathematical model for thin shells filled with liquid. 
The best agreement between theory and experiment was obtained for liquid 
depths greater than the radius of the shell, indicating that surface dynamics 
should be included for lower liquid levels. 

It should not be overlooked that this study attempts to explain only one 
of the physical phenomena associated with the vibration of liquid filled 
elastic tanks subjected to axial excitations. In fact, as demonstrated in 
the experimental phase, it is often difficult to isolate this particular phe- 
nomenon. This analysis is useful, however, in problems where the primary con- 
cern is the force transmitted across the tank bottom (i.e., longitudinal 
dynamics of a missile) and not the local motion of the fluid flow (sloshing) . 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Nov. 3, 1 966 
124-11-05-03-21 
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Figure 2.- Frequency parameter. 
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Figure 8.- Nickel tank natural frequency versus water level. 
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Figure 9 .- First symmetric mode frequency versus water level. 
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Figure 11.- Fundamental symmetric mode frequency response, nickel tank. 
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